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renectance increases inui.uiumcm.y w,m u F uu. <* ..—5 — — modeled after clouds which are uniform in both the vertical and horizontal, retrievals ot 

unity, though not necessarily exact unity for bidirectional reflectance. real-world clouds can be considered the plane-parallel, homogeneous equivalent 

The useful absorbing bands are centered around 1.6, 2.2, and 3.7 pm, with the quantities. That is, the retrievals give the parameters of a homogeneous cloud having 

longer wavelength bands having the greater absorption. These will be referred to nearly the same reflectances as the measured cloud for the particular geometry. Various 

collectively as near-infrared bands (also referred to as shortwave and midwave infrared implementations of this basic approach with the 1.6 and 2.2 pm bands are described in 


» S’ a 

s I z 

s I. 5 

3 O 00 
8 S3 * 


•*,8 1 

o s g 

5 i S 

— 2 o 


C5 <~> 

3 £L 

c „ 


rt' a 
B £ a 

5 p * 

a » U 

0 CTQ O 

1 | & 

| a s 
I -5 1 


d 5. g 


I 3' C/3 

8. « o 

n> a* D5 


H 8- s 

s i §■ 

s & & 

O £D 
O £3 & 

Q- OQ 

o p o 


3' ^ 

> c?. 

•“ I era o 


S g*. S 


■5 s s- 

3 i s 

a* 5 w 

CD _ t/1 


s u i 

| 3 * 
! 5? 

O =• S 

§• & §. 
C -O Si 
Cl EX ^ 


3 o 

o EL s- 


e 2 3 

§ S *o 

8 i g 


^ 00 ' 


O 3 


K S 3 

5 Cl ^ 

CL 0 a 

§ 9 « 

3 2- 3 


o V> CD 

°* 5*. g 

-1 ((a « 


5C ^ < 

2 p o 
3* 9- - 


I* § “ 

I a § 

£ 63 

b: & 2. 

o ^ e 

=r o ^ 

S * S 

n> o cl 


c; fp o 

o* o a* 

2 c £ 


£ § 
cr o p 
a -a 5 
3* g cl 

X3 — O 


may lead to errors or biases in the retrieved cloud effective radius. Clearly, there are other Due to varying absorption, microphysical retrievals made with different near- 

potentially significant and more fundamental sources of retrieval error, such as infrared bands will sample different vertical portions of the cloud and inter different 

uncertainties in the calculation of homogeneous cloud reflectances, instrument error, and effective radii. A retrieval in any given band represents the average radiative effect of the 

atmospheric conditions, to name a few [e.g., Han et al , 1995; Han et al. , 1994; Platnick droplet size profile for that band. A reflectance measurement in the 1.6pm band is least 

and Valero , 1995]. There is also an unexplained tendency for size retrievals to be affected by absorption and can contain information regarding droplet sizes in relatively 


-I § 2 

g. o- 3 


1 5 ~ 


Jb 3" S' 


e x £ 

» o S. 

O ^ = 

g ? & 

Q. O M 


o. n 
o y o' 

Cl & 


3 o 

C «-* 
rt p. -j 


TJ ^ S 

p — ^3 
cs o 
<L % X 


p p 
^ n> =p 
cx n> 


g qo 'f 3 
?T ~ ^ 

^ o < 

£ o ry 

3 •§ 2. 
a e» s- 
s' 8. * 


h 3 (K)' 

o *o a- 

o p o 

2o- 3 

O' 3 oo 


3 P © 
0^3 


(IQ O fl 

OQ C 53 

« era 2 
55. =r 3 


p p ^ 
n* <5 S* 


o ^ 

£ 8 o 
cf 2- oo 

S ^ 

o o 'll 


the literature [Nakajima and King, 1990; Rawlins and Foot , 1990; Twomey and Cocks, significantly larger than in situ measurements, i.e., anomalous absorption [Nakajima et 

1989]. The 3.7 |im band retrievals are complicated by cloud thermal emission which can al., 1991; Twomey and Cocks, 1989], It is within this overall context that effects of cloud 

be a potentially significant part of the measured radiation. Since cloud emissivity will inhomogeneities should also be examined. A better understanding of retrieval uncertainty 

vary' with droplet absorption, thermal emission is a function of droplet size in addition to is especially important given the interest in global droplet size climatologies [Han et al., 

temperature. An algorithm must therefore search for the effective radius that gives a total 1994; King et al., 1992], 
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penetrate only as deep as the z c + dr layer. Photons penetrating deeper would be absorbed 
by the black surface, while photons which have all scatterings above the z c level do not 
contribute to dR. Therefore, the ratio dR/R(r c +dz) must represent the fraction of all 
r e (T)w x {T t T c )dr . (3) reflected photons that penetrate to a maximum depth between z c and r f +dr. This must 

also hold true for arbitrary layers within a cloud. That is, (R{z+ dz)-R(z)) /R{z) gives the 
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4, Multiple scattering weighting functions 


ratio of intensity (radiance) at the top of the cloud scattered into the M ,±A/i/2 bin direction reflectance-based retrieval. In a Monte Carlo calculation, this weighting can be found 

due to incident intensity from the p,± Ap/2 direction (p representing the cosine of the from the total number of scatterings encountered by reflected photons in each layer. The 

zenith angle). Typically, 10 /r-bins of width 0. 10 are sufficient for capturing the angular normalization is then the total number of reflectance scatterings in the cloud. Dividing the 

distribution of intensity in the multiple scattering problems of interest. Bidirectional number of scatterings in each layer by the total number of reflected photons gives an 

reflectance is proportional to S.Ju, Note that n is equated with the cosine of the solar alternative expression for the weighting. Now the normalization becomes the average 
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must be proportional to the number of reflected photons having a scattering in the The accuracy of the proposed weightings in estimating reflectance-based retrieved 

infinitesimal layer. The ratio dfl(T,r f )//?(r r ) therefore gives the fraction of reflected effective radii from vertically inhomogeneous clouds, via Eq. 3, was tested on the cloud 

photons having a single scattering in the infinitesimal layer, a number much less than models described in Table 2 for cloud effective radii increasing with height. Two 

one. As discussed above, this is also equivalent to the average number of scatterings in specifications for the vertical size profile, profiles B (adiabatic) and D (r*~ r) from fable 

the layer for all reflected photons. Integrating this ratio from r to r+Ar then gives the 
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had minor imp^t, typically modifying reeved effective radii by less than 0.2 pm for the are within 0.3 pm for the 1.6 pm band, and within 0.1 pm for the 2.2 and 3.7 pm bands. 

1 .6 and 2.2 pm bands; no differences were found in 3.7 pm size retrievals because of the These differences are small compared with estimates of size uncertainty based on 

band’s greater absorption. The 3.7pm band results use the reflectance signal only, i.e., homogeneous cloud retrievals [ Platmck and Valero , 1995]. The sub-adiabatic aloft 

cloud emission in the band can be considered to have been removed without error. profile results in slightly smaller differences (not shown). Profile D provides the most 

Further, all retrievals are made in the absence of an atmosphere and with a black surface demanding test of the weighting function. Now, size differences are withm 0.9, 0.2, and 
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8 1 8 4 and 9 6 pm, respectively. The estimates for the two longer wavelength bands are single level measurements in validating retrievals from a vertically inhomogeneous 

found to be within 0.2 pm of the retrievals (same analysis as Table 3), whereas the 1.6 pm cloud. Of course, this discrepancy could increase or decrease significantly depending on 

size estimate is about 0.8 pm larger. As in Table 3, the 1 .6 pm band weighting estimate the cloud optical thickness and droplet size profile. Droplet size measurements made at 

tends to differ most from the retrieval. This is due to the band’s relatively weak droplet other cloud levels would differ from the 1 .6 and 2.2 pm band retrievals as well. As the 

viewing and/or solar zenith angle increases, the cloud level consistent w ith the retrievals 
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most, 0.04 pm larger than retrievals made with the black surface, for all bands and ^ ^ ^ gjven by , he first moment of Eq . 5 , which is as 

profiles. expected. Note that r < 0.5 (maximum for overhead solar and viewing angles). This close 

These examples suggest that the effect of a known surface reflectance on droplet proximity to cloud top implies that single scattering retrieval methods will infer larger 

size retrievals from vertically inhomogeneous clouds is negligible. If generally true, droplet sizes than the reflectance-inferred retrievals for adiabatic clouds, except in the 
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surface will differ s.gnificantly from retnevals for the same cloud overlying a black calculations of w. for a cloud overly, ng a black surface can still be used in est.mating 

surface. The assumption of the surface reflectance being known ,s very important. Error retrieved droplet sizes for clouds overlying vegetation, 

due to imperfect knowledge of the surface reflectance is not being investigated here, but 
rather the impact of the surface on the vertical weighting and size retrieval. It is expected 

that reflected photons having a scattering with the surface will increase the information 5.4 Retrievals Based on Single Scattering 
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band. Retrievals from glory or polarization single-scattering reflectance patterns infer 
droplet sizes from the very uppermost region of the cloud (within meters of cloud top), 
and thus may be substantially different from total reflectance retrievals (dominated by 
multiple scattering). For the same reason, single scattering retrievals would be difficult, if 
not impossible, to validate. In addition to providing information regarding the scale of 
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retrievals, or else the retrieval is understood to be the plane-parallel, homogeneous vertical transport, the weightings provide a means for investigating and understanding the 

equivalent value If clouds are inhomogeneous, then separate near-infrared band possibility of an inversion for the droplet size profile using the three near-intrared size 

retrievals may infer different effective radii, a situation which has been observed in retrievals. Analysis shows that the most likely possibility for realizing three unique pieces 

airborne radiometer data [Platnick el at , 1999], We have relaxed the vertical of information from the retrievals is through the use of multiple viewing angles, 

homogeneity constraint in this paper and looked at the effect of modeled vertical droplet Though the emphasis of this work has been on multiple scattering weighting 
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channel. Calculated for a cloud with a total optical thickness of 8, effective radius varying 
as from 5 pm at cloud base to 12 pm at cloud top with profile C in Table 1, cosine of 
solar zenith and viewing angles of /io= 0.65 and /i=0.85, respectively, and an azimuthal 
average. Though the curves are significantly different, the weighting-derived effective 
radii (Eq. 3) differ by only a few tenths of a micron for the profile chosen. 
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Fig. 1 Example of the four analytic models for the vertical profile of effective radius, r e , 
considered in Table 1 for the same prescribed boundary conditions (total cloud optical 
thickness, r c , is 8 and effective radius is 5 pm and 12 pm at cloud base and top. 


Table 1. Analytic models for the vertical structure of effective radius, r e , and liquid 
water content, LWC, versus geometric cloud height, z, and optical depth, r . The following 
convention is used: z increases with height with z- 0 at cloud base to z-h at cloud top, 
optical depth increases towards cloud base with r=0 at cloud top; r c is the total cloud 
optical thickness. The constants ao, aj are found from the optical thickness and the 
prescribed boundary conditions for droplet size at cloud top and base, r, op and r bas e, 
respectively. Other constants (b, c) can be found in a similar manner. 
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Table 3b. Same as Table 3a, but for plane-parallel cloud layers with effective radii decreasing towards cloud top. 
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Table 4. Weighting-derived effective radius retrieval versus cosine of the viewing zenith angle, p. Calculated for r f 8, effective radii 
varying from rt, ase = 5 pm to r lop — 12 jim with an adiabatic profile, and Po _ 0.65. 
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Table 6b. Same as Table 6a, but for plane-parallel cloud layers with effective radii decreasing towards cloud top, and an emission 
weighting, w e , proportional to dr/g*. 
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Fig. 1 . Example of the four analytic models for the vertical profile of effective radius, r e , considered in 
Table 1 for the same prescribed boundary conditions (total cloud optical thickness, T c , is 8 and effective 
radius is 5 um and 12 Jim at cloud base and top, respectively). The top plot shows effective radius as a 
function of optical depth, r. The bottom two plots show the corresponding profile of r t and cloud liquid 
water content, LWC, as a function of geometric height z, where h is the total thickness. The constraints used 
in deriving the profiles are indicated along side each plot. Note that for an otherwise identical cloud 
processes, cloud top r e and LWC would be smaller for the sub-adiabatic profile (shown for *=0.75). 
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Fig. 2 . Two proposed normalized vertical weighting functions, w m (proportional to 
maximum photon penetration) and % (proportional to number of photon scatterings), 
for a 2.2 (im channel, using both superposition formulae (lines) and Monte Carlo 
calculations (symbols). Calculated for a cloud with a total optical thickness of 8, effective 
radius varying from 5 (im at cloud base to 12 pm at cloud top with profile C in Table 1, 
cosine of solar zenith and viewing angles of ^io= 0.65 and fi= 0.85, respectively, and an 
azimuthal average. Plots of the scattering-based weighting function include weightings 
for both reflected and transmitted photons, w' v and w.v, respectively. 
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vertical weighting for transmittance 


Fig. 3. The dependence of the normalized vertical weighting functions on cloud optical 
thickness, z c , for reflectance(top plot) and transmittance (bottom), for a 2.2 pm channel. 
For effective radius varying from 5 pm at cloud base to 12 pm at cloud top with profile C 
in Table 1, cosine of solar zenith and viewing angles of fio = 0.65 and fi= 0.85, 
respectively, and an azimuthal average. 
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Fig. 4. Normalized vertical weighting w m for bidirectional reflectance, for visible and 
near-infrared cloud remote sensing channels, calculated using the superposition formulae 
discussed in the text, and the cloud described in Fig. 2. The cloud optical depth 
corresponding to the retrieved radii for each near-infrared channel is also indicated. 
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Fig. 5a. The dependence of the normalized bidirectional reflectance weighting w m on the 
cosine of the viewing angle, fi, for the three near-infrared channels, and the cloud 
described in Fig. 2. The weighting for reflected flux, or albedo, is also shown. 
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Fig. 6. Example of normalized reflectance and emissivity weightings for a 3.7 (im 
channel. Calculated for a cloud with a total optical thickness of 8, effective radius varying 
as from 5 (im at cloud base to 12|am at cloud top with profile C in Table 1, cosine of 
solar zenith and viewing angles of /io=0.65 and 11 = 0.85, respectively, and an azimuthal 
average. Though the curves are significantly different, the weighting-derived effective 
radii (Eq. 3) differ by only a few tenths of a micron for the profile chosen. 



